clc
clear all

t=0.001:0.001:1
omega_a = 2*pi*1;
omega_b = 2*pi*1;
A = 1;
B = 1;
for i=1:1:1000
    for j = 1:3
        
    end
    vel(i) = A*sin(omega_a*t(i))+B*cos(omega_b*t(i));
       
end

figure (1)
plot(t,vel)


